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ABSTRACT 

We have determined the QCD corrections to the production of squark-antisquark pairs in 
pp coUisions at the Tevatron. If the next-to-leading order corrections are taken into account, 
the renormahzation/factorization scale dependence of the theoretical prediction for the cross 
section is reduced considerably. The higher order corrections increase the production cross 
section at the Tevatron by about a factor two if we compare the next-to-leading order 
prediction at a scale near the squark mass with the lowest order prediction for which, in 
the experimental analyses, the scale was identified with the invariant energy of the parton 
subprocess. This results in a rise of the experimental lower bound on the squark mass from 
the Tevatron by about 20 GeV. 



The search for Higgs bosons and supersymmetric particles rank, sub specie aeternitatis, 
among the most important experimental endeavors of high energy physics [Q. Both areas are 
closely related through the hierarchy problem. In fact, the idea of a supersymmetric exten- 
sion of the standard model is qualitatively supported, though not proved, by the successful 
prediction of the electroweak mixing angle. The colored particles, i.e. squarks and gluinos, 
of supersymmetric theories can be searched for most efficiently at the Tevatron pp collider 
and, in the future, at the LHC up to very high mass values. Lower bounds on the squark 
and gluino masses have been set by both the Tevatron experiments CDF and DO. At the 
90% CL, the masses are larger than about 148 GeV for gluino masses below 400 GeV[|, |]. 

The evaluation of the experimental data has been based so far on the lowest-order (LO) 
production cross sections 0, |^, 0], Figs. |l|(a) and (e). In order to improve the theoretical 
predictions for the pp cross sections, the QCD corrections have to be evaluated in next-to- 
leading order (NLO). The theoretical analysis is straightforward but very tedious due to the 
large variety of production channels and final states. We will report in this letter, as a first 
indispensable step, on the production of squark- antisquark pairs: 

pp^qq + X for q^i (1) 

The qi / qn splitting, particularly important for stop particles, as well as the analysis of gluino 
final states [^] will be discussed in a more comprehensive report at a later time. For the sake 
of simplicity we have taken all squark states mass degenerate, i.e. the rig = 5 squark states 
in the final state as well as the stop particles, which appear in internal loops. The only free 
parameters are therefore the masses of the squarks and gluinos, rrig and rrig, respectively. 



The top mass is fixed to = 174 GeV |10]. 



The next-to-leading order QCD calculation has been performed in the Feynman gauge 
and the singularities have been isolated by means of dimensional regularization. The masses 
have been renormalized in the on- shell scheme. The massive particles are decoupled smoothly 
for momenta smaller than their masses within the modified MS scheme |lll|. When remov- 



ing the infrared divergences, a cut-off A has been introduced for the invariant mass of the 
squark-gluon system in the final state, which separates soft from hard gluon radiation 0. 
If both contributions are added, any A dependence disappears from the total cross section 
for A — >■ 0. The remaining collinear mass singularities can be absorbed into the renormal- 



ization of the parton densities |T2[ , carried out in the MS factorization scheme. The GRV 



parametrizations of the parton densities [13] have been adopted, which allow for proper LO 



and NLO comparisons; CTEQ2 |]14| and MRSH [|T5] parametrizations have also been used 



for assessing the uncertainties from the parton densities. 

If the squarks are lighter than the gluinos, squarks can also be decay products of on-shell 
gluinos, g ^ qq etc.. Since we focus in the present analysis on the evaluation of the QCD 
corrections, we shall restrict ourselves to irreducible final states in which squarks do not 
evolve from on-shell gluinos; after gluino final states are included explicitly, this technical 
assumption will dissolve eo ipso. For the wedge rriq > rrig we disregard, in the same sense as 
above, the decay of the squarks to gluinos. 

Since quark-quark pair initial states give negligible contributions to the generation of 
high invariant mass states in pp collisions, the cross section is built up primarily by quark- 
antiquark annihilation and, to a lesser extent, by gluon fusion. 
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Quark-antiquark initial states. The subprocesses qq —>■ qq are described to leading 
order, Fig. |l](ai) and (02), by mixtures of amplitudes appropriate to unlike/like flavors and 
helicities. Typical standard QCD and supersymmetric vertex corrections are displayed in 
Figs. |l|(&i) and (62)- Ordinary gluon radiation and a diagram related to the renormalization 
of the parton densities are exemplified in Figs. |I|(c) and (d). Note that qq pairs can be 
generated in qg collisions only in next-to-leading order. 

The diagrams have been evaluated analytically to obtain the double differential cross 
sections ddij/dtdu at the parton level; j are the parton indices q, q and t,u are the usual 
Mandelstam momentum transfer variables. The total cross sections atj may be expressed in 
terms of scaling functions fij, 



a. 



nix 



(2) 



They depend on the invariant parton energy through rj = s/Am'^ — 1, and on the ratios 
of the particle masses r = mg/rriq, rt = nit /nig. as is the QCD coupling constant. Renor- 
malization and factorization scales are identified for the sake of simplicity, fiR = fip = Q- 
fij^ denotes the lowest-order contributions. While f^g\ the scaling function associated with 
qg collisions, is zero for qq production, a compact expression |^ can be found for quark- 
antiquark initial states. 



108(1 + r/) 




48r2(l + 77) 



+ Anq(3' + 2(3 



l + r2)2 + 4/5V2(l + r7) 



+ 2/3 + 



(3{l + r,) 



(3) 



108(1 + r/) 



1 + 



+ 2(3 



QLr - 24 + 



48r2(l + 77) 



1 + r2)2 + 4/5V2(l + 77) 



(4) 



where Lj. = log(x+/x_) with x± = 1 + r'^ ± 2(3/(1 ^ (3) and (3 



Amllsfl'^. The 



scaling functions f^^ and , describing the next-to- leading order corrections, are displayed 
in Fig. H. The scaling functions f'^ are split into the "virtual + soft" part (V+S) and 
the "hard" part (H) into which the infrared log-* A (j = 1,2) singularities of the (V+S) 
contribution are absorbed so that these functions are insensitive to the choice of A ^ 0. 
Near the threshold, which is the kinematically most important region in collisions, the 
Sommerfeld rescattering contribution due to the exchange of Coulomb gluons between the 
slowly moving qq pairs, leads to a singularity ~ T^a.sl{3 which neutralizes the phase space 
suppression near threshold. 



with 




(5) 



9(l + r2)2 
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The same relations hold for unlike flavors near threshold. The log^ (3 terms, generated by ini- 
tial state gluon radiation near the threshold, can be exponentiated |]16|. The plateaus for large 
parton energies are due to flavor-excitation and gluon-splitting mechanisms [cf. Fig. [^(ci)]. 
The t- and u-channel exchanges of gluons lead to an asymptotically constant cross section, 
the scale of which is set by the squark mass, i.e. a ~ a^/m?, to be contrasted with the scaling 
behavior a ~ a^/s of the cross section to lowest order. The values of the scaling functions 
in the asymptotic plateau region can be calculated analytically, f^gH ~^ 2159/(194407r) and 
^ -ll/(3247r). 



Gluon-gluon initial states. To lowest-order the diagrams contributing to the subprocess 
99 ~^ s-re the well-known diagrams from quark-pair production [e.g. Fig. 0(ei)] supple- 
mented by the seagull term Fig. |l|(e2) for scalar squarks. Carrying out the next-to-leading 
order program, we find, in the same notation as above 0, 



The gg scaling functions are displayed in Fig. |^(d) in leading and next-to-leading order. Near 
the threshold, the Sommerfeld enhancement is observed again in next-to-leading order. 



with 



(0)t/ir/ \ ^ 7nq7rf3 

99 w; • 



The high-energy plateaus are generated as before by flavor-excitation and gluon-splitting 
mechanisms, /^JJ^ 2159/(43207r) and /j^) ^ -ll/(727r). 

While for ~ 0.1 the higher-order corrections are suppressed significantly with respect 
to the Born term for quark initial states, the large color charge of gluons leads to corrections 
of the order of the Born term for the gg initiated subprocess. 

After these introductory remarks we present our final results in Fig. ^ and Table 0. The 
cross sections of the various subprocesses have been convoluted with the parton densities in 
the GRV |]13[, CTEQ2 and MRSH parametrizations. The following conclusions can 



be drawn from the analysis of the figures and the table. 

(i) It is obvious from Fig. 3(a) that the theoretical predictions for the pp qq production 
process are improved considerably by taking into account the next-to-leading order QCD 
corrections. While the dependence on the renormalization/factorization scale Q is quite 
steep and monotonic in leading order, the Q dependence is significantly reduced in next-to- 
leading order for reasonable variations of the scale, running even through a broad maximum 
near Q ~ m^/3. Since the cross section is built up mainly by the quark channels [^85%] 
and thus based on well-measured parton densities, the variation between different parton 
parametrizations is negligibly small. 

(ii) The K factors, defined as the ratio K = (Tnlo/o'lo [with all quantities in the numer- 
ator and denominator calculated consistently in NLO and LO, respectively], depend only 
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mildly on the squark and gluino masses. Experimental mass bounds can therefore be cor- 
rected easily for higher-order QCD effects. A sample of K factors is collected in Table |l| for 
various choices of the renormalization/factorization scale parameter Q. 

(iii) Finally in Fig. 3(b) we illustrate the impact of the QCD corrections on the exper- 
imental lower bounds of the squark masses. We compare the lowest-order cross section at 
Q = \/S, the scale adopted in experimental analyses, with the next-to-leading order pre- 
diction at Q = rriq [^. The NLO cross section is significantly larger at this theoretically 
reasonable scale than the LO cross section at the scale a/^. Taken at face value, this increases 
the bound on the squark mass by about 20 GeV. While no precise value of Q can be defined 
a priori, it is clear nevertheless that the experimental bounds derived from Tevatron data 
0, ^ are very conservative and that the true bounds are likely to be higher by as much as 
~ 20 GeV. 

We thank S. Lammel for useful discussions on the Tevatron squark and gluino mass limits 
and their theoretical basis. 
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Table 1: K factors for a set of g, g masses and a range of renormalization/factorization scales 



Q at the Tevatron energy ^/s = 1.8 TeV; GRV parton densities [jl3| 



Figures 



Figure 1: Generic diagrams for squark-antisquark pair production: (oi), (02) basic Born-level 
qq diagrams; (61) and (62) vertex corrections due to gluon and gluino exchange; (c) gluon 
emission; (d) gluon-quark process; (ci), (62) gluon fusion. 



Figure 2: The scaling functions for squark-antisquark pair production in (a) qq, (b) q'q, (c) 
qg and (d) gg collisions. The notation follows eq.(0) with t] = s/irri^ — 1; V+S denotes the 
sum of the virtual and soft corrections, H the contribution of hard gluon emission. Mass 
parameters: niq = 250 GeV and nig = 200 GeV. 



Figure 3: Total cross section for the irreducible production of squark-antisquark pairs 
pp — > qqX at the Tevatron energy ^/s = 1.8 TeV. (a) Dependence on the renormaliza- 
tion/factorization scale Q for the leading order (LO) and the next-to-leading order (NLO) 
predictions, and sensitivity to different parton densities; mass parameters as in Fig. 2; (b) 
Dependence of the cross section on the squark mass for rrig = 200 GeV; GRV parton den- 
sities for NLO, and EHLQ parton densities for LO used in the experimental analysis. [^. 
Upper full line of the NLO prediction corresponds to the renormalization/factorization scale 
Q/rriq = 1/3, middle line = 1 and lower line = 2. 
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